Quantum Dynamics of Radical-Ion-Pair Reactions 
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Radical-ion-pair reactions were recently shown to represent a rich biophysical laboratory for the 
application of quantum measurement theory methods and concepts, casting doubt on the validity of 
the theoretical treatment of these reactions and the results thereof that has been at the core of spin 
chemistry for several decades now. The ensued scientific debate, although exciting, is plagued with 
several misconceptions. We will here provide a comprehensive treatment of the quantum dynamics 
of radical-ion-pair reactions, generalizing our recent work and elaborating on the analogy with the 
double-slit experiment having partial "which-path" information. This analogy directly leads to the 
general treatment of radical-ion pair reactions covering the whole range between the two extremes, 
that of perfect singlet-triplet coherence and that of complete incoherence. 



I. INTRODUCTION 

Spin-selective radical-ion-pair reactions are a rare ex- 
ample in chemistry where spin degrees of freedom and 
their relatively small interaction energy can have a dis- 
proportionately large effect on the outcome of chemical 
reactions. The study of radical-ion-pair reactions is at 
the core of spin chemistry [l| , by now a mature research 
field directly related to photochemistry Q and photosyn- 
thesis 

Radical-ion pairs are biomolecular ions created by a 
charge transfer from a photo-excited D* A donor-acceptor 
molecular dyad DA, schematically described by the reac- 
tion DA -> D*A -> D* + A*~, where the two dots rep- 
resent the two unpaired electrons. The magnetic nuclei 
of the donor and acceptor molecules couple to the two 
electrons via the hyperfine interaction, leading to singlet- 
triplet mixing, i.e. a coherent oscillation of the spin state 
of the electrons. The reaction is terminated by the re- 
verse charge transfer, resulting to the charge recombi- 
nation of the radical-ion-pair and the formation of the 
neutral reaction products. It is angular momentum con- 
servation at this step that empowers the molecule's spin 
degrees of freedom to determine the reaction's fate: only 
singlet state radical-ion pairs can recombine to reform the 
neutral DA molecules, whereas triplet radical-ion pairs 
recombine to a different metastable triplet neutral prod- 
uct. If the lifetime of radical- ion pairs is large enough, 
it is seen that hyperfine and Zeeman interactions, neg- 
ligible compared to thermal energy, can produce a large 
effect on the reaction yield. 

Theoretically, the fate of radical-ion-pair reactions is 
accounted for by the time evolution of p, the density ma- 
trix describing the spin state of the molecule's two elec- 
trons and magnetic nuclei. The master equation satisfied 
by p currently is at the center of a scientific debate since 
it was shown [i[ that from the advent of spin chemistry 
until now, the master equation traditionally used to pur- 
sue all theoretical work was phenomenological, masking 
the existence of non-trivial quantum effects and leading 
to an incomplete understanding of radical-ion-pair reac- 
tions. A fundamental master equation was derived based 
on quantum measurement theory, as the radical-ion-pair 



recombination process was interpreted Q to be a con- 
tinuous quantum measurement of the radical-pair's spin 
state. 

Since almost all possible theoretical predictions that 
can be made on radical-ion-pair reactions are founded 
on the master equation for the radical pair's density ma- 
trix, it is clearly important for the scientific community 
to converge to a common understanding of the quantum 
foundations of spin selective radical-ion pair reactions. 
Such reactions determine the late-stage dynamics in pho- 
tosynthetic reaction centers 5, 6], and furthermore there 
is increasing evidence that radical-ion-pair reactions un- 
derlie the avian compass mechanism, i.e. the biochemical 
compass used by migratory birds to navigate through the 
geomagnetic field |7H13j. A deeper understanding of the 
quantum dynamics inherent in radical-ion-pair reactions 
will thus enhance and quite probably significantly change 
our understanding of these biological processes. 

Recent works have, however, challenged the theoretical 
description developed in [4|, sparking a debate on which 
is the fundamental master equation describing the spin- 
state evolution of reacting radical-ion pairs. In order of 
appearance, Jones & Hore [l4j have produced a different 
master equation for p, conjectured to follow from quan- 
tum measurement theory. Shushin [l5j claimed that the 
Bloch-Redfield relaxation theory already accounts for the 
quantum dynamics of radical-ion-pair reactions and no 
new concept from quantum measurement theory is ac- 
tually needed. Most recently, Ivanov et. al. [16| and 
Purtov [l7| put forward a derivation of the radical-ion- 
pair master equation, leading to the traditional equation 
of spin chemistry. 

We will here address in detail the misconceptions con- 
tained in these recent papers and then move on to ex- 
amine in depth the quantum measurement dynamics of 
radical-ion pair reactions, attempting to shed light on 
several aspects of this system and their physical inter- 
pretation. We will also address our own misconceptions, 
showing that our previous work Q describes part of the 
picture, and a more general theoretical description is nec- 
essary. Such a description has been developed in [l8[ 
where it was shown that radical-ion pair reactions are 
analogous to the archetypal quantum paradigm, the dou- 



2 



D*A 




a 



Open Quantum System 

Radical-Ion-Pair 

— — — — — — "\ 

magnetic i 
(*D + *A") S interactions (*D + "A") J k T 



Triplet 

Measurement 
Apparatus 



(DA) T 



FIG. 1: Energy level structure and reaction dynamics of 
radical-ion pairs. A photon excites the singlet neutral precur- 
sor molecule DA into D*A, and a charge transfer creates the 
radical- ion pair. The excited vibrational levels (DA)* of the 
neutral DA molecule form the measurement reservoir, which 
has two functions: it acts as a measurement device for the 
spin state of the radical-ion pair, and it acts as a sink of 
radical-ion pairs, i.e. if a radical-ion pair recombines into the 
singlet channel, the electron tunnels into a reservoir state , 
and a fast spontaneous decay results in the ground state DA 
(which is the singlet product) and a photon emission. Similar 
for the triplet reservoir. 



ble slit experiment with partial " which-path" informa- 
tion. In the following section we outline the physical 
system under study, while in Section III we analyze the 
conceptual problems of the approaches in [l4j], [H| 
and [17|. In Section IV we present a detailed analysis of 
the general theory, and among other things, we analyze 
single molecule experiments in order to elucidate severe 
inconsistencies of the traditional and the Jones- Hore the- 
ories. Finally in Section V we address the issue of energy 
conservation in radical-ion-pair reactions, which serves 
as another fundamental benchmark in the comparison of 
the various theoretical descriptions. 



II. RADICAL-ION PAIRS AS AN OPEN 
QUANTUM SYSTEM 

We will first describe in detail the physical system un- 
der consideration. The energy levels and reaction dy- 
namics of radical-ion pairs are depicted in Fig[TJ We 
do not consider the dynamics of creation of radical-ion 
pairs and we assume a given initial condition for the 
pair's spin state. For example, in many cases radical- 
ion pairs are created in the electronic singlet state and 
the nuclear spins arc completely random, i.e. they are 
in the maximally mixed state. What we care about in 
this work is the combined dynamics of (i) the coher- 
ent singlet-triplet mixing induced by the magnetic in- 
teractions and (ii) the irreversible spin-state-dependent 
charge recombination that terminates the reaction. For 
the study of the quantum dynamics of the recombination 



process, we consider a radical-ion-pair with its center of 
mass fixed. In other words, we completely neglect ef- 
fects stemming from diffusion in solutions, collisions with 
other molecules, spin relaxation etc, since at this point we 
care to understand the fundamental quantum dynamics 
inherent in the radical-ion pair. Now, in all experimen- 
tal studies we have a macroscopic number of radical-ion 
pairs. Each one of them is a single open quantum system 
regarding the spin degrees of freedom. It is an open quan- 
tum system because the recombination process inherent 
in each molecule disturbs a what would be a unitary spin 
evolution. Moreover, due to recombination, radical-ion 
pairs disappear from the ensemble and are replaced by 
neutral chemical products. That is, what we deal with 
is an ensemble of open quantum systems with varying 
number of particles. This non-trivial aspect of radical- 
ion pair reactions and the confusion stemming thereof 
is in part fuelling the recent debate on the fundamental 
master equation describing these reactions. 



III. THEORETICAL DESCRIPTION OF 
RADICAL-ION-PAIR REACTIONS 

Since the advent of the field of spin chemistry, the spin 
chemistry community has described the evolution of the 
spin state of radical-ion pairs by the master equation 
(from now on to be called " traditional" ) 



Qsp)~^t(pQt + Qtp) (1) 



where the unitary evolution of the radical-pair's spins 
state is driven by the magnetic interactions embodied in 
H, and the other terms attempt to represent the reaction 
taking place through the singlet (rate kg) or the triplet 
(rate kx) recombination channel. A few definitions are 
in order. The dimension of the density matrix and all 
operators acting on the radical-ion pair's Hilbert space 
is An, where 4 is the spin multiplicity of the two elec- 
trons and n = (2Ii + l)(2i2 + l)-..(2/fc + 1) is the nuclear 
spin multiplicity of the k magnetic nuclei (residing in the 
donor and acceptor molecules) and having nuclear spins 
Ii, J2,...,Ifc. The operators Qs and Qt are the singlet 
and triplet projection operators, given by 



Qs 
Qt 



1/4 

3/4- 



si ■ s 2 
si ■ s 2 



(2) 



These two projector operators add up to unity, Qs + 
Qt = 1, i.e. the radical-ion pair is either singlet or 
triplet. 

In 2009 Kominis 0] argued that this master equation 
misses a lot of the physics of the recombination process. 
By properly identifying the quantum system and the en- 
vironment "watching" the system's evolution, the master 
equation describing the evolution of the radical-ion-pair's 
spin state until it reacts was derived in [4| : 



d Pn 



dt 



-i[n, P ] 



ks + k T 



(pQs + Qsp-2QspQs) (3) 
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This is a trace preserving master equation that describes 
the decoherence of the radical-pair's spin state brought 
about by the internal to the molecule recombination dy- 
namics (the subscript in dp nr denotes the "not reacted" 
molecules). This equation applies to a single quantum 
system. However, we have to model reactions, i.e. dis- 
appearance of radical-ion pairs and creation of neutral 
products. That was modelled by another equation evolv- 
ing the existing number of radical-ion pairs, 

dN = ~Ndt(k s (Q s ) + k T (Q T )) (4) 

which is interpreted as follows: in the time interval be- 
tween t and t+dt the probability of the radical-ion pair to 
be in e.g. the singlet state is (Qs) — Tr{pQg}, and hence 
the probability of recombination is ksdt(Qs) ■ Since there 
exist N{t) radical-ion pairs, N(t)ksdt(Qs) of them will 
disappear into singlet products, and similar considera- 
tions apply to the triplet channel, leading to the update 
rule for the number N(t) shown in As will be ex- 
plained in the following, this description is part of the pic- 
ture, and applies to states with maximal singlet-triplet 
coherence. The general case of incoherent mixtures is 
more complicated. In 2010 Jones & Hore [14J introduced 
a different master equation derived from a rather ad-hoc 
application of the theory of generalized quantum mea- 
surements and positive maps. The equation reads 

dp/dt = -i[H,p] - k s (pQ s + Q SP - QspQs) (5) 
- k T (pQ T + Q T p - QtpQt) 



with some kind of apparatus. The measurement is inter- 
nal in each molecule, as has already been explained [4]. 
When a number of radical-ion pairs have reacted, the re- 
action products and the number thereof represent purely 
classical information. Whether or not we acquire this 
classical information is completely irrelevant and cannot 
have any back-reaction whatsoever on any quantum sys- 
tem. Finally, the authors have not derived their mas- 
ter equation (obviously because it is not derivable) us- 
ing the usual recipe of considering the specific Hamilto- 
nian interactions between system and environment and 
then tracing out the environment. Their ad- hoc postu- 
lated master equation following from a hand-waving use 
of the operator-product formalism completely lacks phys- 
ical foundations. Moreover, as will be shown in Section 
V, the Jones-Hore theory is plagued by a severe inconsis- 
tency: if one applies it to a single- molecule experiment, 
the average of all possible realizations does not reproduce 
the predictions of the Jones-Hore master equation. The 
traditional approach faces the same problem, and both 
theories also violate conservation of energy, as will be 
shown in Section VI. 

Most recently another three contributions appeared in 
the literature. Ivanov et. al. claimed to have "consis- 
tently" derived the traditional master equation, and Pur- 
tov has solved exactly a particular example and shown 
that the solution is what would follow from the tradi- 
tional master equation. What both works have in com- 
mon is the starting point of their derivation, which rep- 
resents a deep physical misconception. 



A. Conceptual problems of the Jones-Hore Theory 

The conceptual problem with the Jones-Hore theory is 
the erroneous application of single-quantum-system mea- 
surement dynamics to a many-particle ensemble of inde- 
pendently reacting molecules. In particular, the philoso- 
phy of the Jones-Hore derivation is that in a box with 
a macroscopic number of radical-ion pairs, the detec- 
tion of the particular reaction products during the time 
interval dt affects the quantum state of the remaining 
molecules. Specifically, these authors argue like this: if 
ksdt molecules react to produce singlet products, this will 
have the effect of projecting the rest of the system onto 
"a more triple" state, i.e. they density matrix acquires 
a term ksdtQxpQr- This is an unphysical statement, 
for the following reasons: first, the molecules in the box 
are completely independent and non-interacting systems. 
The fact that molecule i reacted cannot in any way affect 
the quantum state of molecule j. Second, in the quantum 
measurement of single quantum systems, one usually con- 
siders the quantum system and the measurement appara- 
tus (another quantum system) coupled to it. Due to this 
coupling, the two become entangled, and of course any 
information extracted from the apparatus back-reacts on 
the system. In the case of radical-ion-pair reactions, we 
do not perform a global measurement on all molecules 



B. The derivations of Ivanov et. al. and Purtov 

The starting point of both derivations is the asser- 
tion that the radical-ion pair quantum state, call it 
\D + A~) and the neutral product state, call it \DA), 
are eigenstates of one and the same Hamiltonian Ho, 
and transitions are induced between them by some 
perturbation V. That proposition, namely that reac- 
tants can be in a quantum superposition with the re- 
action products of an exothermic irreversible chemical 
reaction is utterly wrong. One can think of a simi- 
lar situation, the nuclear decay of uranium into bar- 
ium and krypton, and attempt to write for the total 
state — cu|U) + CB a |Ba) + CK r |Kr), and continuing 
along these lines one could arrive at states of the form 
| nuclear bomb exploded) + 1 nuclear bomb not exploded). 
Or, to give an example familiar to chemists, one could 
argue along the lines of Ivanov et. al. that molecular 
oxygen, molecular hydrogen and water are eigenstates 
of one and the same Hamiltonian, and a perturbing po- 
tential causes transitions between them. If that would 
be the case, Nature would be dramatically different and 
we would not be reading this text now. Thus, the con- 
ceptual problem with these two derivations is that the 
authors build up a model that does not correspond to 
a physical entity. Hence their result is of no practical 



4 




FIG. 2: Radical-ion-pair reactions pictured with Feynman di- 
agrams. (a)The reaction as modelled by Ivanov et. al. and 
Purtov (b) The reactions as modelled by Kominis (c) Analogy 
with nuclear fission. It is the dense manifold of intermediate 
states that is ignored in the "effective" theory of Ivanov et. 
al. and Purtov. Similarly, it is the intermediate state of 236 U 
that prohibits the formation of any quantum coherence be- 
tween the parent nucleus and the fission product. 



of physical reality. Philosophically, this viewpoint treats 
the whole universe as a single Hilbert space and induces 
coherent transitions between disjoint parts of this space. 
If nature indeed behaved like this, Schrodinger cats would 
abound around us, in fact we ourselves would not exist in 
our current form to discuss this. We remind the reader 
that in the cold atom community various groups have 
created coherent oscillations between molecules and free 
atoms (see for example [Hj]), but at nano-Kelvin temper- 
atures needed for Bose-Einstein condensation. We are far 
(if ever) from doing the same in electron transfer or, for 
that matter, nuclear reactions. 



C. Shushin's Argument 

Shushin [l5| has argued that the Bloch-Rcdficld re- 
laxation theory already explains radical-ion pair reaction 
dynamics and no new concept is required. The Bloch- 
Redfield theory is a second-order perturbation theory, 
and indeed similar terms like QsP or QspQs appear there 
as well. However, the theory describes decoherence in- 
duced by a noisy environment. This has nothing to do 
with the charge recombination process. This is a funda- 
mental decoherence process internal in the molecules. We 
cannot suppress this decoherence, no matter how good 
(i.e. noise free) experiments we perform. Shushin points 
out the main aspect of our previous work, namely the 
spin decoherence of not-yet-reacted radical-ion pairs, but 
regards it as an ad hoc interpretation. We claim that 
this interpretation directly follows from the formalism 
developed in Q, and its physical origin will be further 
elaborated upon in the following Section. 



IV. TOWARDS THE COMPLETE THEORY 



interest. Essentially, the authors build a so-called "ef- 
fective" theory, ignoring a dense manifold of interme- 
diate states and the instant decoherence caused by the 
decay of these states. In Feynman diagram terminology, 
these authors look at the radical-ion-pair reaction as a 
single vertex joining \D + A~) with \DA) and a phonon, 
as shown in Figj2£i. The actual radical-ion-pair reaction, 
however, that takes place in nature, is described by the 
process shown in Fig|2j), where the radical-ion pair is 
converted into a vibrational excited state of the neutral 
molecule, (DA)*, which then emits a photon and we then 
reach the ground state DA, which is the reaction prod- 
uct. The photon instantly "expels" to the environment 
the information on it's precursor state, and hence any co- 
herence between the reactants and the reactions products 
can not be physically sustained and hence it is pointless 
and erroneous to consider it theoretically. The analogy 
with nuclear fission is shown in FigJ5£. In other words, 
to pertain that there exists a quantum state of the form 
\t/j) = a|D + A~) + /3\DA) is a completely distorted view 



In [l8[ we have outlined the derivation of the complete 
quantum dynamics description of radical-ion-pair reac- 
tions, while elucidating an analogy between this system 
and the double slit familiar from quantum optics. We 
will here elaborate on the complete theory in more de- 
tail. One part of the complete theory is the consistent 
treatment of the reaction. The other is the treatment of 
the not-yet-reacted radical-ion pairs, which undergo spin 
decoherence. 



A. Decoherence of Nonreacted Molecules 

What the master equation ([3]) says is that the internal 
measurement going on in each radical-ion pair is respon- 
sible for spin decoherence. This internal measurement 
measures the observable Qs with a total measurement 
rate k$ + kr, as explained in A process contributing 
to loss of spin coherence can be nicely visualized by a 
Feynman diagram like the one shown FigEJ That is, a 
virtual tunneling from the radical-ion pair state to the 
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FIG. 3: A virtual process contributing to the decoherence of 
the radical-ion-pair spin density matrix. This decoherence is 
described by the density matrix equation derived by Kominis 
Q. In this virtual process, driven by the interaction Hamilto- 
nian in the second order perturbation theory, the radical-ion 
pair couples to the excited vibrational state, which, instead of 
decaying to the ground state, couples back to the radical-ion 
pair. These processes are unobservable and the source of spin 
decoherence of nonreacted molecules. 

singlet reservoir state and back has the effect to dephase 
the coherent-singlet triplet oscillation. Similar is the ef- 
fect of the triplet reservoir. Such virtual processes are 
unobservable. This is the physical interpretation of trac- 
ing out the reservoir degrees of freedom, while assuming 
a still-existing radical-ion pair, as explained in [4|. 



state |S')(S'| because we found a singlet product? Not 
necessarily. To claim that is a gross misunderstanding of 
the state reduction. Consider for simplicity a single two 
level atom observed by a unit efficiency photo-detector. 
If the detector clicks, i.e. a photon is observed, we know 
for sure that the atom is (now, just after the click) in 
the ground state \g). We know nothing as to what state 
the atom was in before observing the photon. That is, 
we cannot say that just before observing a photon the 
atom was for sure in the excited state |e). The atom 
could have been at any quantum superposition of the 
form a\g) + (3\e). In other words, the observation of 
the photon updates our knowledge about the state of the 
atom, it does not inform us about the atomic state prior 
to the observation. Similarly, the observation of a singlet 
product molecule does not at all mean that just prior 
to the observation the radical-ion pair was in the singlet 
state. It might have been in any coherent superposition 
of singlet and triplet states, and the molecule decided to 
recombinc through the singlet channel. 

Suppose now that all N radical-ion pairs are at time 
t prepared in the same state \ip) — a\S) + fi\T), hence 
the single-molecule density matrix is p\ = |a| 2 |S')(S'| + 
\f3\ 2 \T)(T\ +aj3*\S) (T\ + a*(3\T)(S\. The ensemble den- 
sity matrix is p — Npi. Clearly, in the time interval be- 
tween t and t + dt wc will detect dns — ksdtTr{pQs} — 
k s dtN|a| 2 singlet and dn T = k T dtTr{pQ T } = k T dtN|/3| 2 
triplet neutral products. The change in p will then be 
dp = pidN + Ndp\. That is, p changes because radical- 
ion pairs react (first term) and because non-reacted 
radical-ion pairs are subject to a state change. From 
the work of Kominis Q it follows that this state change 
is singlet-triplet decoherence, and it is described by ([3]). 
Let us now see what each of the theories has to say about 
dp. 



1. Kominis Theory 



B. Maximum Coherence Extreme: Problem with 
the Traditional and Jones-Hore Theory 

In all subsequent discussions we neglect nuclear spins 
and focus on a two-level radical-ion pair, the Hilbert 
space of which is spanned by two states, the singlet 
and the triplet \T) (we should say the triplet |To), the 
zero projection state of the triplet manifold spanned by 
\T±) and \Tq), but for notational simplicity we write \T)). 

Suppose wc have a macroscopic number N of radical- 
ion pairs. Each is described by a density matrix pi, hence 
the density matrix of the ensemble is p ~ Yli=i Pi- Sup- 
pose that particle j recombines into the singlet channel. 
What is the effect of this recombination on Rl Clearly, 
now R = Pi- That is, we have to remove one par- 

ticle from the ensemble of particles, and from the total 
quantum state R we have to remove pj. What was the 
state of the particle to be removed? Was it the singlet 



From (U]) we find that dN = —dn$ — dnx, whereas 
from ([3]) we find that dpi = —(ks + kT)dtp co i 1 /2, where 
p coh = N(a/3*\S)(T\+a*f3\T)(S\) is the coherence of p. 
Hence in this theory we find 

dp = -(dn s + dn T )pi - — ^-^-^ jO co h (6) 

2. Traditional Theory 

From |T]) it is straightforward to find that 

, /j j \ 1 ,dn s + dn T (k s + k T )dt 
dp = -(dn s + dn T )pi + ( — )p 

- \{(dns - dn' s )\S)(S\ + (dn T - dn' T )\T)(T\) (7) 
where dn s = /ctM 2 an d dn T = ks\(3\ 2 . 
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3. Jones-Hore Theory 
From ([5]) we find that 
dp = -(dn s + dn T )pi + ^ dn S +^ dn T _ ^ + p 
+ dn s \S)(S\+dn T \T)(T\ (8) 



channel, for example we take ks — 0. What we naturally 
expect to find is that half of the radical-pairs will disap- 
pear into triplet products, leaving us at the end with a 
density matrix p^ = ^-\S)(S\ describing the other half 
which have remained in the non-reactive singlet state. 

1. Kominis Theory 



4- Physical Problem of the Traditional and the Jones-Hore 
Theories 

The common problem of the traditional and the Jones- 
Hore theory is in the second term of their expressions 
for dp. Namely both theories update the quantum state 
of the nonreacted molecules based on the recombination 
events of those molecules that reacted. This unphysical 
state of affairs has been analyzed in [l8[ in analogy with 
the photon double slit experiment. We briefly reiterate 
the main point: the fact that in the time interval dt we 
have dns singlet and driT triplet recombinations affects 
the density matrix of the non-reacted molecules, as if 
there is a weird communication between radical-ion pairs 
in the ensemble, is unphysical, as unphysical is to assert 
that the detection of the number of photons on the ob- 
servation screen in Young's double slit experiment affects 
the quantum state of the other photons flying through the 
slits. In an ensemble of 1000 say radical-ion pairs being in 
a coherent singlet-triplet superposition, the fact that 10 
of them recombined into the singlet channel in the time 
interval between t and t+dt cannot have any influence on 
the spin state of the other 990 whatsoever (unless they 
are entangled, which is a situation far from possible at 
the current experiments). In contrast, according to the 
theory of Kominis, the nonreacted molecules suffer loss 
of their spin coherence as time goes on, which is a sin- 
gle molecule effect acting internally and independently in 
each molecule. 



C. Maximum Incoherence Extreme: Problem with 
the Kominis Theory 

A maximally coherent state is one extreme possibility. 
The other extreme is the maximally incoherent state. We 
will study this case now and then move on to the more 
realistic general case. Suppose that at some instant in 
time we know that the ensemble density matrix is an 
equal mixture between singlet and triplet, i.e. suppose 
that 

Pt = ^\S)(S\ + ^\T)(T\=N Pl (9) 

where Pl = (\S)(S\ + \T)(T\)/2. That is, according to 
Pt there are N/2 radical-ion pairs in the singlet subspace 
and another N/2 in the triplet. We again consider the 
absence of magnetic interactions, so there is no singlet- 
triplet mixing, and we consider just one recombination 



As before, dp — pidN + Ndp\. Since p\ = (Qs + 
Qt)/2 it easily follows that dp nl — i.e. there is no more 
coherence to be lost from an incoherent mixture, hence 
the master equation Q for nonreacted molecules does 
not induce any change in p. Since furthermore (Qs) = 
(Qt) = 1/2, the equation determining the number of 
radical-ion pairs (Q} gives dN = —kTdtN/2, hence we 
get dp/dt = —krp/2 i.e. p^ — (meaning all radical-ion 
pairs disappeared into products), which obviously is not 
the case. So it is clear that the machinery of removing 
particles according to (0} does not work in this case [20| . 
The reason is that, writing p t — Np, upon measuring dnx 
triplet recombination products we should update p t by 
the rule pt+dt — Pt + dp, where dp = —dn^Qri whereas 
Kominis' theory updates pt by dp — —dn^p. In other 
words, in this extreme incoherent case we know for sure 
that the observation within the time interval dt of dnr 
triplet products implies that these products came from 
triplet state radical-pairs, so our knowledge about the 
ensemble state is updated to a state having "less triplet 
character" , i.e. we have to remove from p the projection 
dnrQrpQT- 

2. Traditional and Jones-Hore Theories 

It is straightforward to show that both the phenomeno- 
logical as well as the Jones-Hore theories gives 

dp = -k T dt^^- (10) 

leading to p^ — NQs/2, as it should be. Thus in this, 
maximum incoherence extreme, both of these theories 
work fine. 

V. GENERAL CASE 

In [l8[ we developed the formalism addressing the gen- 
eral case, which is this: suppose that at time t the radical- 
ion-pair ensemble is described by the density matrix p t , 
which can be any mixture. In the following time interval 
dt we will surely find dns — ksdtTi{pQs} singlet and 
dnT = /cTrfiTr{pQx} triplet neutral products. How can 
we consistent with all information at hand update pt into 
Pt+dt = Pt + dp? For completeness, we reiterate the an- 
swer to this question. If p t were an incoherent mixture of 
singlet and triplet radical-ion pairs, then we could inde- 
pendently remove "singlet" and "triplet" character from 
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p according to dpi nco h = -ksdtQ s pQs - k T dtQ T pQT- 
In the other extreme, if p t was a state with maximal 
singlet-triplet coherence, then we would have to remove 
the complete information on the state of the reacted 
molecules, which is what we did in our previous work 
0: dp C oh = -(dns + dn T )pt/Tr{pt}, where p t /Tr{p t } is 
the single-molecule density matrix. In the general case 
we have at our disposal a measure of singlet-triplet co- 
herence of p t . This measure was defined in [lj| as 



Pcoh = 



Tr{/3 S TPTs} 
Tr-tpsSj-TMpTT} 



(11) 



dp 
~dt 



dp„ 



where p S s = QspQs, Ptt = QtpQt, Pst = QspQt 
and pts = QtpQs- The incoherent part of p is pincoh = 
Pss+Ptt, while the coherence is p co h = PST+Pts- More 
on p co h in the following. Finally, we also have to add 
the change in the density matrix of the not-yet-reacted 
radical- ion pairs, dp nl -, given by ([3]). Putting everything 
together we get 



At 



(1 - Pcoh){ksQspQs + k T QTpQr) - Pcoh(k s Tr{Qsp} + k T Tr{Q T p}) 



Tr{p} 



(12) 




tQcf 1 ) 



FIG. 4: Radical-ion-pair recombination dynamics in the ab- 
sence of singlet-triplet mixing and ks ~ 0, starting from the 
initial state + \T))/s/2. (a) Time evolution of Tr{Q s p} 
and (b) of Tr{Qxp} as predicted by (|12[) (red solid line) and 
the traditional as well as the Jones-Hore theory (blue dashed 
line) . 



A. Reaction of Spin Coherent Radical-Ion Pairs 

To illustrate the fundamental differences between the 
various theories, we used [l8j a very simple scenario of 
an ensemble of radical-ion pairs all starting out in the co- 
herent state (\S) + |T))/\/2, with no singlet-triplet mix- 
ing (i.e. we zero the magnetic Hamiltonian) and just 
one recombination channel (e.g. we set ks = 0). It is 
straightforward to find that both the traditional as well 
as the Jones-Hore theory predict that half of the rad- 
ical pairs will react (of course producing triplet prod- 



ucts) while the other half will remain locked in the non- 
reacting singlet state (Fig|4|. In glaring contrast, the 
newly developed theory predicts that 75% of the rad- 
ical pairs will react, and 25% will remain in the sin- 
glet state forever. At first look this sounds completely 
counter-intuitive and almost unacceptable, according to 
the layman understanding of quantum mechanics. Since 
the initial state is a coherent superposition of singlet 
and triplet with equal amplitudes and the singlet is non- 
reacting, how on earth do 75% of the molecules react? 
The resolution of this apparent paradox is this: the tra- 
ditional theory is based on a derivation like the one pre- 
sented by Ivanov et. al. [jfjj . This theoretical construct 
allows the following evolution: Starting with a radical- 
ion-pair state -^{\S) + |T)), and because the singlet is 
non-reactive and the triplet reactive, this state evolves 
to -^( | singlet radical pair) + | triplet product)) and of 

course that results in 50% triplet products!! But such an 
evolution is impossible. As we explained previously, the 
reaction does not create coherent superpositions between 
reactants and reaction products. Asking an ill-defined 
question in quantum physics can lead to paradoxical re- 
sults. Simple projective-measurement quantum mechan- 
ics says that when we measure a physical observable A in 
a state which is a coherent superposition of the observ- 
able eigenstates, say \tp) = ci|ai) + C2I&2) we obtain a\ 
with probability |ci| 2 and so on. So a question that can 
be asked regarding the initial state -^(\S) + \T)) is this: 
what is the probability in a projective measurement of 
Qs to find zero or one? Of course both probabilities are 
1/2. But this does not tell us anything about the reaction 
products. This is so because the radical-ion pair might be 
in this coherent state at time t and there is a finite prob- 
ability of reaction within the following time interval dt. 
But the radical- ion pair might choose not to react, and 
its spin state will then evolve according to ©, hence the 
final product yield is not that straightforward to predict 
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by using the naive projective measurement argument. To 
reiterate, there is no observable having a radical-ion-pair 
state and a neutral product state as eigenstates. This is 
the source of the confusion and the apparently paradox- 
ical result we mentioned before. 



B. Single Molecule Scenarios 

We will here analyze single molecule experiments to 
further illuminate the deficiencies of the traditional and 
Jones-Hore theories. Consider a single radical-ion pair at 
time t in the coherent superposition state 775 (1-5) + \ T)). 
Again, we eliminate singlet-triplet mixing and take k$ = 

0. Clearly, (Q s )t = (Qr)t = 1/2. Within the next 
time interval dt, there is a finite probability dpT — 
kTdt(Qx)t = k^dtjl of recombination in the triplet chan- 
nel. In a particular realization of this experiment that 
probability will materialize, so at time t + dt we have no 
radical-ion pairs, hence (Qs)t+dt = (Qr)t+dt = 0. How- 
ever, both the traditional as well as the Jones-Hore the- 
ory predict (Qs)t = 1/2 for all times. How are these two 
contradictory facts reconciled? Well, we have to consider 
a different realization of this single molecule experiment, 
in which the radical-ion pair does not react within the 
time interval dt after its preparation in the coherent state 
at time t. How does the state of an nonreacted molecule 
evolve according to the Jones-Hore theory? In their 
derivation of the Jones-Hore master equation [14|, the 
authors argue: during the time interval dt we will have 
ksdt radical pairs recombine through the singlet channel, 
making the rest more "triplet", i.e. adding to p the term 
kgdtQTpQr- Similarly we will have another k^dt radical 
pairs recombine through the triplet channel, adding to p 
the term kxdtQspQs- And finally, we add to p the term 
(1 — ksdt — kxdt)p because a fraction 1 — ksdt — kxdt 
of radical-ion pairs did not recombine at all. Thus 
the update rule for the Jones-Hore master equation is 
Pt+dt = (1 - k s dt - k T dt)p + ksdtQrpQr + k T dtQ s pQs, 
which leads to (0. So according to this theory the den- 
sity matrix of a non-reacting radical-ion pair will change 
by dp = 0. In other words, the state of a non-reacting 
radical- ion pair remains -^(\S) + \T)) until the time it 
reacts, i.e. d(Qs) = 0. But in the previous realization 
we had d(Qs) = —1/2, so to keep d(Qs) = 0, we must 
have realizations with d(Qs) > 0, so that the average 
satisfies the prediction of §5§ shown in FigQJi. So it is 
clear that the Jones-Hore theory is not self-consistent, 

1. e. it cannot account for single molecule scenarios in 
such a way that the average of them reproduces the pre- 
dictions of the Jones-Hore master equation ([5]). To re- 
solve this problem and keep d(Qs) = the Jones-Hore 
theory performs a rather unphysical act, which necessi- 
tates the existence of more than one radical-ion pairs: 
in an ensemble of molecules N, the fact that a number 
of them reacts through the triplet channel projects the 
rest to a state with more singlet character. Although we 
are repeating ourselves here, it is a crucial point and we 



will reiterate: the fact that in our box containing 10 23 
radical-ion pairs a number of them (say 10 6 ) recombined 
within the time interval dt to create triplet products, this 
very fact projects the rest 10 23 — 10 6 molecules in a state 
that is more singlet than the state they were at time t. 
We claim that this is an impossible act. To reiterate, the 
Jones-Hore theory fails to account for single-molecule re- 
alizations in a self-consistent way, and manages to keep 
d(Qs) = by some weird action at a distance between 
independent molecules in an ensemble. The traditional 
theory does not even in principle consider what happens 
to a non-reacting molecule. Such a concept does not exist 
in the traditional theory, therefore the above considera- 
tions apply to the traditional theory as well. 

In contrast, the picture emerging from our description 
is the following: Take kxdt = 1, so that starting with the 
coherent state -^(\S) + \T)) the probability for triplet 
recombination is kxdt/2 = 1/2. We will trace the state 
evolution at time steps dt, as shown in Fig(5] In a par- 
ticular realization of this single molecule experiment this 
recombination takes place, and the reaction terminates. 
In another realization the radical-ion pair will not react in 
this time interval dt after its preparation in the state 1^), 
so according to ^ spin coherence will be lost. At time 
t = we had p co h = 1, but if the radical- ion pair does not 
react, it will lose half of its coherence (easily found from 
©) and at time t — dt the coherence measure p co h will 
be reduced to p co h — 1/4. At this point there is again 
the possibility of triplet recombination, the probability of 
which is still 1/2. Again, if this recombination does not 
take place, coherence will be still reduced, and at time 
t = 2dt we will have p co h = 1/16. All these are possi- 
ble realizations, as is the the realization that we have no 
reaction until the coherence has been reduced to exactly 
zero, where we have an incoherent single-triplet mixture. 
At this point the radical-ion pair will either remain in the 
non-reacting singlet state (probability 1/2) or recombine 
in the triplet channel (probability 1/2). So the longest 
the reaction can proceed is until dp nl /dt — 0. The ques- 
tion is, however, when is the coherence really zero? In 
other words, how many realizations do we have to average 
in order to estimate the triplet yield? Again p co h is the 
measure of when the reaction is really over. So the triplet 
yield Y T is not p tr + PnrPtr + pl r Ptr + — = 5 + ^ + — + ^r, 
but each term in the series is suppressed by p co h: i-e. 

Y T = Ptt + Pcoh(0)p nr Ptr + Pcoh(^)PnrPtr + ■•■ (13) 

This is how the result shown in FigfJJi is produced. 



C. More on p co h 

The meaning of the definition of p co h is rather ob- 
vious if we consider a mixture of N c radical-ion pairs 
in the coherent state a\S) + /3\T), Ns radical-ion pairs 
in the singlet state \S) and Nt radical-ion pairs in the 
triplet state \T). The density matrix of such a mixture 
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p = y|S><S| +i|T><T| 



S><T| + y|T><S 




Pcoh(0)=l 



triplet 
p,=l/2 ' product 



time 



p = ^SxS| +i|T><T| 
+ ^S><T|+^T><S| 

p = y|S><S| +i|T><T| 
+ i]SxT| +-i|T><S| 

p = i|S><S|+ijT><T| 



P m =l/2 
Pcoh(*)=l/4 



P„r = l /2 

p coh (2dt)=l/16 



P„=l/2 



Pc„h(-)=0 



= 1/2 



triplet 
product 



triplet 
product 



triplet 
product 



=1/2 



singlet 
product 

FIG. 5: Single radical-ion-pair state evolution with no singlet- 
triplet mixing, a single recombination channel (ks = 0) and 
a coherent initial state. The probability for triplet recombi- 
nation is ptr = 1/2 an d the probability for no recombination 
iS Pnr = 1/2. 



is p = (N s + N c \a\ 2 )\S)(S\ + (N T + iV c |0| 2 )|T>(r| + 
N c a/3*\S)(T\ +N c a*f3\T)(S\. In this case Tr{p S TPTs} = 
N 2 | Q | 2 |/3| 2 , Tr{p SS } = N s + N c |a| 2 and Tr{p TT } = 
N T + N c |/3| 2 . As intuitively expected, if iV c = we 
have p co h = 0, and if Ns = = we get p con = 1. 
Now let us take into account the nuclear spin degrees 
of freedom. For simplicity we consider just one nuclear 
spin with spin 1/2 and we denote by |±) the nuclear 
spin basis states. If the radical-ion-pair spin state is 
for example -^(\S) + |T))|+), then again p con = 1. If, 
on the other hand the nuclear spin is in an incoher- 
ent mixture of |+) and |— ), then the spin state of the 
radical-ion pair will be described by the density matrix 
p=l(\S)(S\ + \T)(T\ + \S)(T\ + \T)(S\)(\+)(+\ + \-)(-\) 
and we find that p co h = 1/2. In this case the electronic 
spin is in a maximally coherent state for which we ex- 
pect p co h = 1 but the nuclear spin is in mixed state, and 
therefore p co h = 1/2 underestimates the electron spin co- 
herence. This is because p co h depends on TrjpsTPTs}, 
which is also a measure of the mixed character of the 
spin state. However, it is just electron spin coherence 
that we care about regarding the recombination process. 
We therefore have to normalize TrjpgxPTs} by Tr{p 2 }, 
and accordingly we normalize Tr{pss} and Tt{ptt} by 
Tr{p}. Now p co h = 1 for the aforementioned mixed nu- 



clear spin state, as it should be. 



VI. ENERGY CONSERVATION IN 
RADICAL-ION-PAIR REACTIONS 

In radical-ion-pair reactions, singlet-triplet mixing 
conserves energy, that's why the magnetic Hamiltonian 
mixes the singlet with the degenerate To state of the 
triplet manifold. It is well known that the exchange inter- 
action, of the form — Jsi • S2 splits the singlet and triplet 
manifolds by J and hence suppresses singlet-triplet mix- 
ing. We study this effect and compare how the three 
theories deal with energy conservation. We consider a 
radical-ion pair with a single spin-1/2 nucleus and take 
as an initial state | -M — ), i.e. the two electrons (first two 
entries) are in the T + state and the nuclear spin I is in the 
|— ) state. We consider an isotropic hyperfine coupling 
asi ■ I of one electron with the single nucleus and we also 
consider the exchange coupling between the electrons, 
hence the magnetic Hamiltonian \sT-L = — Js\ ■ S2 + asi ■ I. 
Since during the reaction we lose radical-ion pairs into re- 
action products, we consider the quantity 



AE(t) = Tr{ Pt H} - Tr{p t=0 H}Tr{p} t 



(14) 



That is, we consider the expectation value of the sys- 
tem's energy, Tr{pH}t, which obviously tends to zero 
due to the reaction removing molecules, minus the en- 
ergy the system would have at time t if the energy per 
molecule remained the same. Stated otherwise, AE(t) is 
the change in the energy per radical-ion-pair times the 
number of radical-ion pairs. Since in our considerations 
regarding the fundamental dynamics the radical-ion pair 
does not exchange energy with the environment, we ex- 
pect AE(t) = 0, at least within 1/t, i.e. energy should be 
conserved, or at least, if it is not, this non-conservation 
should be within the limits allowed by Heisenberg's un- 
certainty. The difference AE(t) is given units of the ex- 
change coupling J, whereas time has units of 1/J. The 
product AE(t)t should thus not exceed unity. However, 
this is the case with the traditional and Jones-Hore theo- 
ries, in contrast with the result following from the general 
theory developed here and in [18| . This is shown in Fig|6l 
It is seen that the Jones-Hore theory (the same result 
holds for the traditional theory) violates energy conser- 
vation by 5<7, whereas our description keeps the system's 
energy change within the physical limits. 



A. Analogy with Quantum Dots 

A similar issue comes up in the system of two cou- 
pled quantum dots interrogated by a point contact, suc- 
cinctly analyzed by Gurvitz (2lj . The tunneling rate in 
the point contact, T, depends on whether the left or the 
right dot is occupied, while a coherent oscillation takes 
place between the two extreme states (a) left well occu- 
pied and (b) right well occupied. Using Bloch equations, 
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t (104/ J) 

FIG. 6: Evolution of AE(t) predicted by (a) the Jones-Hore 
theory and (b) this work, i.e. equation (|12l) . We consider 
a radical- ion pair starting in the state | + H — ), i.e. the two 
electrons are in sth spin up state and the single nuclear spin in 
the spin down state. We take ks = 0, kr = 1.0, the magnetic 
field is zero, the hyperfine coupling with the single nucleus is 
a = 1.0 and the exchange coupling J = 20.0. 

Gurvitz describes the decoherence brought about by the 
measurement device of the double-well coherent oscilla- 
tions. There are three rates that determine the problem: 
the coherent oscillation frequency Qo, the energy mis- 
match e of the ground state energies of the two wells 
and the decoherence rate T, derived from the point con- 
tact tunneling rates. When the decoherence rate T = 0, 
the coherent oscillations are evident in the probability 
of occupation of the left well. For increasing T, the co- 
herent oscillations are suppressed and the density ma- 
trix tends to the fully mixed state (in this case it is a 
two-dimensional density matrix and hence a aa —> 1/2 
as t —> oo). The first point of relevance to radical-ion 
pairs is that, as Gurvitz notes, this is a counter-intuitive 
manifestation of the quantum Zeno effect, namely that 
the effect of a strong continuous measurement is to de- 
localize the state of the measured system, as opposed 
to the more familiar localization produced by a series of 
projective measurements. The exact analog with radical- 
ion pairs is the decoherence of the nonreacted molecules. 
The evolution dp m /dt brings about exactly such a de- 



localization. The second point is that there is a sub- 
stantial probability to occupy the right well (50% as 
t — > oo), even though the energies of the two wells dif- 
fer by e 7^ 0. That is, the electron tunneling from the 
left to the right well seems not to conserve energy. This 
apparent paradox is easily resolved, since the measure- 
ment apparatus, the quantum point contact, can readily 
exchange energy with the quantum system. Clearly so, 
because the interaction Hamiltonian coupling the left well 
with the point contact is essentially an electrostatic inter- 
action, embodied in the "Coulomb-blockade" -like term 
"Hint =~ SflirC^CLiajar + h.c.) where c L (c£), a; [a]) 
and a r (a\) are annihilation (creation) operators for the 
left well, the left reservoir and the right reservoir of the 
point contact barrier, respectively. In other words, the 
electron on the left well interacts with the point con- 
tact and can hence make a transition with AE 7^ 0. In 
radical- ion pairs, however, the coupling with the "envi- 
ronment" is the tunneling analyzed in 0], the result of 
which is a pure decoherence. This is so because there 
is no mechanism for the radical-ion pair to exchange en- 
ergy with the vibrational reservoir, hence the spin state 
evolution must conserve energy. 



VII. CONCLUSIONS 

We have analyzed in detail two extreme cases, that 
of maximal singlet triplet coherence, and that of maxi- 
mal incoherence. The fundamental theories of radical-ion 
pair reactions that existed so far were applicable only at 
one of the two extremes. We have developed the general 
theory that spans the continuous range between these 
two extremes. Furthermore, we have elucidated the fun- 
damental conceptual problems faced by the traditional 
master equation of spin chemistry as well as the more re- 
cent variant by Jones & Hore. We have unraveled the in- 
consistencies of these theories by attempting to describe 
single-molecule measurement scenarios. Finally we have 
compared all theories in the front of energy conserva- 
tion. In all cases the general theory we have developed 
consistently describes several gedanken-experiments. It 
is clear, as noted in [l8j], that when the singlet-triplet co- 
herence (quantified by p co h) is small, the traditional the- 
ory will provide a satisfying description of spin-selective 
radical-ion pair reactions. Such a case would be experi- 
ments dominated by spin relaxation. On the other hand, 
in experiments suppressing relaxation effects to the ex- 
tent that p C oh is large, we expect significant deviations 
from the previous theoretical understanding. 
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